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We study the eflFect of cooling on the spatial distribution of the topological charge density in 
quenched SU(2) lattice gauge theory with overlap fermions. We demonstrate that as the gauge 
field configurations are cooled, the Hausdorff dimension of regions where the topological charge is 
localized gradually changes from d = 2-^3 towards the total space dimension. Therefore, the cooling 
procedure destroys some of the essential properties of the topological charge distribution. 
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Introduction 

Topological charge density is an important characteris- 
tic of the QCD vacuum, recently involved in phenomeno- 
logical studies of many new hypothetical effects [H-Q]. 
However, the spatial structure of the topological density 
distribution seems to be not well defined since the rele- 
vant properties of the underlying vacuum structure de- 
pend on the measuring procedure 0, IBl- The classical 
instanton approach Q assumes that the nonperturbative 
physics is governed by the scale of Aqcd, which means 
that the dimensionful quantities like volumes occupied 
by topological fermion modes should depend on Aqcd 
but not on the lattice spacing. On the contrary, the lat- 
tice measurements demonstrate that these volumes do 
depend on the spacing (i.e. on the measurement resolu- 
tion) and shrink to zero in the continuum limit [lol - [l^ . 

It turns out that the continuum definition of the topo- 
logical charge density 




FIG. 1: Isosurfaces of the topological charge density 
q{x) = ±10"'* for a fixed time slice, corresponding to the 16 
lattice in Tabled Colors represent positive (red) and negative 
(blue) values, respectively. For the animation, see [l9l ]. 
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cannot be directly applied to the lattice gauge theory, 
since the discretized version of ^ is no longer a full 
derivative. There are two widely used methods to study 
the topology of gauge fields on the lattice. First, one can 
apply a smearing procedure, which makes the gauge fields 
smoother and thus closer to the classical fields. Second, 
one can rely on the lattice version of the Atyah-Singer 
theorem and define the total topological charge of a gauge 
field configuration as the number of zero modes of the 
overlap Dirac operator 15[ on this configuration. The 
corresponding local density of t opologi cal charge can be 
defined, for example, as follows |16l - [l8| : 
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where D(x,x) is the zero-mass Neuberger operator and 
the trace is taken over spinor and color indices. Another 
attractive property of this definition is that it allows us 
to measure a local imbalance in the number of left- and 



right-handed quarks (chirality), which is important for 
lattice studies of the local CP-violation in strong inter- 
actions [1^ . A typical result of the lattice simulation for 
this quantity (without cooling) is shown in Fig. [TJ 

At the moment there are many investigations related 
to the spatial structure of the topological charge distri- 
bution [ll, [13, mHlBl , which use both of the alternative 
definitions. The measurements which rely on the cool- 
ing procedure mostly sug gest an instanton-likc picture 
of the QCD vacuum [2^, while the definition ([2]) typ- 
ically shows that the topological charge is localized at 
low-dimensional objects (defects) [H, [ij, [H, [l^ and has 
a very- long-range structure of the distribution [2^. At 
the qualitative level it is known that both definitions yield 
the top olog ical charge densities which arc strongly cor- 
related |2ll . [27I [28j . For an alternative filtering method 
based on adjoint fermions see Ref. [29| . 

The aim of this paper is to fill the existing gap in the 
literature and to demonstrate in what way the cooling 
procedure affects the dimensionality of regions where the 
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a [fm] 


Li X Lt 


V [fm*] 


# conf 


3.200 


0.117 


12^ X 12 


3.93 


50x13 


3.295 


0.100 


14^ X 14 


3.90 


50x13 


3.332 


0.094 


15^ X 15 


3.89 


50x13 


3.365 


0.088 


16^ X 16 


3.88 


50x13 


3.425 


0.078 


18^ X 18 


3.87 


50x13 



TABLE I: Lattice parameters used in the calculation: cou- 
plings /3, lattice spacings a, lattice sizes x Lt, physical 
volume V , and number of gauge field configurations. 



topological charge density is localized. We use the defi- 
nition ^ based on zero modes of overlap Dirac operator 
and show that as the gauge field configurations are cooled 
the dimension of these regions gradually tends to 4, which 
is the total space dimension. The procedure makes the 
effective resolution of the measurement lower and thus 
provides a result close to the instanton picture. We ver- 
ify our result using several measures of the localization 



iry our resu 



Technical details 



We work in the quenched SU{2) lattice gauge th eory 
with the tadpole-improved Wilson-Symanzik action [3l[ . 
Lattices we used are listed in Table HI We also imple- 
ment the cooling procedure described in Ref. with 
coefficient c = 0.5 for the APE-smcaring. For each lat- 
tice spacing we consider thirteen different stages of the 
cooling procedure: 0, 1, 2, 5 - 12, 20 and 50 iterations 
of the algorithm. For valence quarks we use the Neu- 
berger's overlap Dirac operator [l5|. Its eigenvalues and 
eigenfunctions are given by the following relation 



(3) 



The quantities we measure in the present work are 
functions of two basic ingredients: the "chiral conden- 
sate" computed on a mode with eigenvalue A, 



(4) 



and "chirality" computed on a mode with eigenvalue A 
[in agreement with the definition 



(5) 



Here we sum over spinor and (omitted) color indices. The 
total values of both chiral condensate and chirality are 
given b y a n infinite sum over all eigenvalues. Lattice 
studies |32l . l33j suggest that the long-distance properties 
of QCD can be treated with a finite cutoff of the fermionic 
spectrum. We hereby restrict our consideration to the 
IR part of the Dirac spectrum consisting of zero modes 
(A = 0) and few low-lying modes (A ^ 0). 



Inverse participation ratio (IPR) for an arbitrary nor- 
malized distribution a{x) is usually defined in the follow- 
ing way 



(6) 



wliere is tlie total number of lattice sites x. From this 
definition one can clearly see that IPR = iV if a{x) is 
localized on a single site and IPR = 1 if a{x) = const, i.e. 
the distribution is unlocalized. In general IPR is equal 
to the inverse fraction of sites occupied by the support of 
a{x). Since this fraction of sites can be thought of as a 
number of four-dimensional lattice liypercubes covering 
tlie support, the Hausdorff dimension d of tliese regions 
can be extracted from the asymptotic behavior of IPR at 
small lattice spacings a 



IPR(a) 
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where c is a constant. It is also useful to mention, that 
in physical units IPR^^ is equal to the part of the total 
volume occupied by the distribution. 

In the following sections we will modify the standard 
definition ^ to adapt it to our particular cases (i.e. un- 
normalizcd or non-normalizablc distributions, etc.). The 
final result will show an equivalence of the chosen defini- 
tions. 



Ordinary IPR for zero modes. 

In this section we compute the inverse participation 
ratio for the fermionic zero modes according to the one 
defined in RcL 



12|: 



IPRo = A^ 



X 



(8) 



A=0 



where the brackets [...]a=o denote an averaging over all 
zero modes and further averaging over all gauge field con- 
figurations. Results are presented in Fig. [2l 

The left-hand figure shows how the localization de- 
pends on the lattice spacing a - the finer the lattice, the 
larger the IPR. This fits very well to the idea of vanish- 
ing total volume occupied by fermionic zero modes in the 
continuum limit a — > (see Ref. Q for a review). Using 
the fit ([7]) we recover the fractal (Hausdorff) dimension 
d of the volume. Results for the fits with fixed numbers 
of cooling steps arc presented in the Table HIl Here, to 
minimize errors, we also prepared an alternative sample 
consisting only of those configurations which do not lose 
all the fermion zero modes during the cooling. We picked 
tlien the values with better (and also sufficient) statistical 
significance. 





steps 
2 steps 
7 steps 
9 steps 
1 steps 
20 steps 
50 steps 



0.075 0.08 0.085 0.09 0.095 0.1 0.105 0.11 0.115 0.12 
Lattice spacing (fm) 



16 
14 
12 
10 
8 
6 
4 
2 




Hi 



10 



Number of cooling steps 

FIG. 4: Chiral IPR for the lowest nonzero modes. First definition, Eq. (|10p . 
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Number of 
cooling steps 


Fractal 
dimension 


Standard 
error 


P-value 





2.84 ± 0.44 


15% 


0.008 


1 


2.66 ± 0.66 


25% 


0.027 


2 


2.49 ± 0.46 


18% 


0.013 


5 


2.17 ± 0.49 


23% 


0.021 


6 


2.75 ± 0.66 


24% 


0.025 


7 


3.17 ± 0.51 


16% 


0.009 


9 


3.71 ± 0.34 


9% 


0.001 


12 


3.88 ± 0.23 


6% 


4 • 10-* 



TABLE II: Fractal dimension of the fermionic zero modes and, 
equivalently, of the topological charge distribution. 
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FIG. 6: Fractal dimensions at various cooling stages. The 
solid line is shown to guide the eye. 



Chiral IPR for low-lying modes. First definition. 



Chiral IPR for zero modes. Second definition. 



In this section we modify the IPR to measure local- 
ization properties of the topological charge distribution. 



The average ehirality 



E 



is zero, therefore we 



J A 



have to use either the absolute value |p^(a:;)| or the square 

[p|(a;)] ^ . Here we stick to the definition from [s^l, which 
in our terms has the following form 



IPRg = N 



X 



(9) 



A=0 



Results are presented in Fig. [3] From the plots we 
conclude that the topological charge distribution behaves 
similar to the zero modes, tending to occupy a vanishing 
volume in the continuum limit. We can also compute the 
chiral IPR for small but nonzero eigenvalues (in our case 
we pick first 7 eigenvalues, A < 200 MeV), 



IPRi^o = N 




(10) 



Chiral IPR for these modes is small (Fig. and thus 
the topological charge distribution at this part of the 
spectrum is delocalizcd. 



Finally we consider a second definition of the chiral 
IPR according to [l0| : 



IPRg = N 



Eko(^)r 

X 

E i^o(^)i 



(11) 



A=0 



where, as before, p^{x) denotes the ehirality on a zero 
mode ([2]). Results are presented in Fig.O As can be seen 
from Figs. [51 [31 and [5l the IPR for the zero modes and 
for the topological charge density on these modes are the 
same up to negligible deviations. Results of the fitting 
procedure coincide for these three cases and are shown in 
Table [ill The coincidence is not accidental, because for 
the zero modes \D,^^\ = and 7^|'!/'o) = i|V'o)- This 
means that on a given mode Pq{x) and Pq(x) are equal 
to each other up to a sign. 



Fractal dimension. Results and conclusions. 

To conclude, we demonstrate that the topological 
charge is localized on low-dimensional fractal structures, 
whose fractal (Hausdorff) dimension depends on the 
number of cooling steps. The obtained dimension is 
about d = 2 3 for a few [n < 6) steps of the cool- 
ing, while it grows to d = 4 with further iterations (see 
Fig. [6]). For a long cooling {n > 20) the result becomes 
insignificant, because the procedure leads to a delocal- 
ization of the distributions as can be seen in Figs. [2H51 
(otherwise IPR remains consistent with a constant within 
error bars). We suppose that it can be caused by the an- 
nihilation of the instanton/anti-instanton pairs. Indeed, 
comparing the mean action evolution (Fig. [7]) with the 



5 





1 
1 


1 1 


12? 

1^4 

16^ : 

-1 n4 

18 


o 
CO 








CO 
c 


0.1 






o 








o 
nj 








Mean 








0.01 







20 40 60 80 100 

Number of cooling steps 



FIG. 7: Mean action at various cooling stages. 

one from e.g. Ref. [s^] we see that the annihilation phase 
in our case could start already from n ^ 20. In Ref. p^ . 
where the same cooling algorithm is used, the annihila- 
tion takes place even at a smaller number of steps. 
The main conclusions of our paper are the following: 

(1) Fermionic zero modes and chirality are localized on 
structures with fractal dimension d = 2 3, which 
is an argument in favor of the vortex/domain- wall 
nature of the localization (sH, . 

(2) A long sequence of iterations of the cooling proce- 
dure provides a result close to the instanton picture, 
i.e. destroys the low-dimensional structure of the 
QCD vacuum. 

Finally, let us briefly mention a possible phenomeno- 



logical consequence of our study. One of the most promis- 
ing effects appearing due to the nontrivial topology of 
the QCD vacuum is the so-called "chiral magnetic effect" 
(CME) [H, which states the generation of an electric cur- 
rent in parallel to an external magnetic field. Topological 
charge density in this case can be understood as an imbal- 
ance in the number of left- and right-handed light quarks 
induced by a nontrivial gluonic background. This effect 
is expected to explain charge asymmetries observed at 
RHIC [111. Some evidences of the CME on the lattice 
as well as numerical estimates for the values of the local 
topological charge were also obtained in Refs. (ssl - li^ . At 
the current level of analytic studies CME is considered 
as an effect on the background of spatially homogeneous 
axial fields [4^, while the lattice simulations predict an 
irregular structure of the would-be axial field (see Fig.[T|). 
This spatial inhomogeneity can be treated within a chi- 
ral superfluid model |6| , where the chirality is carried by 
an effective axion-like field. Knowledge of the nature of 
the topological charge localization can help us to trans- 
late lattice Euclidean properties of the chiralityto the 
language of an effective Minkowski field theory [4J1 . 
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